Harmonic Measure for Percolation and Ising Clusters Including Rare Events 
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We obtain the harmonic measure of the hulls of critical percolation clusters and Ising-model 
Fortuin-Kastelyn clusters using a biased random-walk sampling technique which allows us to measure 
probabilities as small as 10~ 300 . We find the multifractal D(q) spectrum including regions of small 
and negative q. Our results for external hulls agree with Duplantier's theoretical predictions for 
D(q) and his exponent —23/24 for the harmonic measure probability distribution. For the complete 
hull, we find the probability decays with an exponent of —1 for both systems. 

PACS numbers: 05.45 Df, 05.40 Jc, 41.20 Cv 



The harmonic measure is a fundamental property of 
geometric objects, and its determination is of consider- 
able interest. For a given object, it may be denned in 
the following way: consider the object to be a grounded 
conductor with fixed charge of unity. Then the harmonic 
measure, /i, is the distribution of electric field on the hull 
(surface) of the object. Alternatively, we may imagine 
that we allow many random walkers to start far away 
from the object and record where they land. The prob- 
ability density of hitting the hull at a point is /i. This 
quantity is of both theoretical [TJ [2] and practical [3] in- 
terest. In particular, if the shape in question is fractal, 
fi shows interesting scaling properties. 

In this Letter we show how to find fi numerically for 
two fundamental systems that produce fractal clusters 
in two dimensions, namely percolation [1] and Fortuin- 
Kasteleyn (FK) [5] clusters in the Ising model. Our 
method allows us to sample very small probabilities (of 
order 10 -300 ) using random walker simulations. In the 
case of these fractals, /i is of particular theoretical im- 
portance because it is multifractal. Our large dynamic 
range allows us to explore this property fully. 

The harmonic measure /i is non-negative and normal- 
ized on the hull: J dfj, = 1. A partition function [5] can 
be defined by dividing the hull into j boxes of length I, 
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where pj = J d/i over box j. For large fractal clusters Z q 
scales as a power-law, 
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where R is the length scale of the cluster and 
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T(q) = (q-l)D(q) 

is the multifractal scaling exponent; D(q) is called the 
generalized dimension. We recall some special values of 
D(q): -D(O) is the fractal dimension of the support of 
the measure, which describes the region the hull covers. 



Additionally, D(l) — 1 is known from Makarov's theorem 
[TJ. A related quantity of interest is the curve f(a), which 
is the Legendre transform of r(q): 
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The rest of this paper will focus on r(q) and D(q) . We 
are also able to produce f(a) by numerically taking the 
Legendre transform of our measured r(q). 

The exact spectrum of D(q) has been obtained for per- 
colation [5] and the more general Q-state Potts model [5], 
derived from generalized conformal invariance in terms of 
a central charge c, 
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Percolation and FK clusters for the Ising model corre- 
spond to c = and c = 1/2, respectively. 

Eq. (0 was derived in the scaling limit for the acces- 
sible or external hull of the clusters [HI HO] • For a finite 
system the external hull is approximately produced by 
closing all fjords on the complete hull with a neck size 
of order unity. The closing of the fjords reduces the di- 
mension for percolation and the Ising model from 7/4 and 
11/8 for the complete hull to 4/3 and 5/3 for the external 
hull, respectively [8]. The theory derived for percolation 
compares well with previous simulation results [T]. How- 
ever, these simulations were limited to relatively small 
systems, » 10 5 sites and did not probe deep into the 
fractal surface. 

The authors of Ref. [1] used the method mentioned 
in the first paragraph: a large number of walkers were 
released far away from the cluster and allowed to diffuse 
until they were absorbed on the hull. This method is only 
able to measure the growth probabilities to an accuracy 
of ss 1CP 10 , while percolation clusters with 10 5 sites can 
have regions of the hull with probabilities smaller than 
1CP 100 . Although these regions do not contribute to D(q) 
for large q, they have a significant impact on D(q) for 
small and negative q. The algorithm of this paper can 
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FIG. 1: The signpost algorithm. Left: the hull of the cluster 
(line) and the adsorbed random walkers (dots). Right: the 
number of walkers that hit each of the sections of the hull (A- 
F). Probe step, a) 18 random walkers are released and adsorb 
onto the hull. The histogram shows the distribution of hits 
with a dashed line representing the probability threshold, 4. 
The large dots on the histogram show where the two ends of 
the signpost should be placed, b) 18 more walkers are released 
and adsorb onto the hull and the signpost (grey), exactly 9 
walkers on each. The histogram shows the probability on the 
hull, c) The next probe step: 18 walkers are released from 
the locations along the signpost where walkers in the step b) 
landed (open circles). These walkers have half the weight as 
the ones released in parts a) and b). The dashed line and 
circles on the histogram show the new threshold and signpost 
end points. Color online 



measure probabilities down to 10 -300 which completely 
samples lattice systems with ~ 10 4 hull sites. We have 
applied this method on systems as large as ~ 10 6 hull 
sites. In this paper we consider only percolation and 
Ising clusters though our method is quite general and 
can also be applied to off-lattice clusters. 

In order to treat very rare events, we use a biased 
sampling of random walkers, in which we keep track of 
the "lucky" random walkers that penetrate deep into the 
fjords of the hull. The method has three steps: In the 
first iteration, N random walkers with equal weight are 
released from outside the cluster and allowed to diffuse 
until they touch the hull. Each walker has weight l/N. 
The weights of the walkers are temporarily added to the 
probability distribution of the site where they land. This 
step probes the cluster to find where the regions of small 
measure lie. The hull sites that bound regions below a 
threshold probability (say 0.1 to start out) are used as 
the end points of an absorbing line, or signposts. Then, 
the probability added in the first step is removed from the 
distribution and N more random walkers are released in a 
similar fashion as the first step. These walkers can either 
diffuse onto the hull or a signpost line. The weights of the 



FIG. 2: The results of our simulation on the complete hull of 
a percolation cluster, where the harmonic measure goes from 
high to low though the colors red, green, then blue. The scale 
is « 50 orders of magnitude. Sites outside cluster are white, 
and inside, grey. Color online. 

walkers that touch the hull in this step are permanently 
added to the probability distribution. After all walkers 
have been absorbed, the signpost lines are removed. 

In the next iteration, the probe step is repeated with 
TV walkers released from the locations along the signpost 
lines where walkers absorbed in the previous iteration. 
The threshold for small probability is reduced by a con- 
stant factor, e.g., 10. These walkers will have their weight 
reduced by the percentage of the walkers in the last step 
that touched a signpost line. This method is similar to 
methods used in chemical physics to numerically obtain 
reaction coordinate information [IT]. Fig. [2] shows the 
harmonic measure of the external hull of a percolation 
cluster obtained using this method. 

Our simulations are performed on a periodic triangular 
lattice with height h and width w such that h = 100u>, 
so that we obtain clusters that wrap around in width but 
not in height. One ambiguity which must be resolved is 
the definition of random walkers touching the hull. Here, 
we interpret this as the walker hopping onto a hull site. 

The percolation clusters are grown using the Leath al- 
gorithm [T2], with p equal to the site threshold for the 
triangular lattice, p c — 1/2. If a given cluster spans the 
width of the system, the top hull of the cluster is found 
using a simple border walking algorithm related to the 
method of generating percolation cluster perimeters by 
random walks |13j . The list of complete hull sites of the 
cluster is then used in the signpost method to obtain the 
harmonic measure. If the topmost vacant sites bordering 
the cluster are used instead of the occupied sites as the 
adsorbing sites, one obtains the external hull. 

For the FK clusters on the Ising model, bonds are 
placed between adjacent same-spin sites with probabil- 
ity p c = 1 — e~P Jc , where on the triangular lattice 
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FIG. 3: D(q) vs. q for the complete and external hulls of per- 
colation clusters for four different widths compared to theory, 
Eq. Q '2 (black line). Inset: small q behavior of D (q). Color 
online. 



p c = 1 — 1/^/3. We use the Swendsen-Wang method 
P3] to equilibrate the system and simultaneously gener- 
ate the FK clusters. After the system is sufficiently equi- 
librated, we attempt to find a spanning cluster. These 
spanning bond clusters must be converted to site clus- 
ters if they are to be used with our algorithm. We do 
this by making another triangular lattice with half the 
lattice spacing, such that all bonds are centered on even 
sites in the new lattice. Bonds are copied to the new 
lattice at the even sites on which they are centered. Odd 
sites are added to the cluster if two adjacent bonds meet 
at the odd site. Next, the perimeter-walk algorithm is 
used to record the locations of the hull sites; then we use 
the signpost method to find the harmonic measure. As 
in percolation, an external hull can be obtained, though 
for Ising clusters it requires the addition of artificial va- 
cancies to all sites bordering the cluster. 

The signpost method iteratively obtains smaller and 
smaller probabilities by reducing the weight of the ran- 
dom walkers released in each round, in our case by a 
factor of 10. We took the number of walkers, N to de- 
pend on the system width, w. For example, for w = 400 
we use N — 2 ■ 10 6 and for w = 4000, N — 2 ■ 10 7 . The 
signpost method is performed until all probabilities have 
been measured or until the minimum measurable proba- 
bility 10 -300 has been reached. This minimum is close to 
the smallest value that can be stored in a double precision 
floating point number. Significantly smaller values could 
be obtained by storing the log of the probability instead 
of the probability itself, though we have not implemented 
it. 

The locations of the sites and their associated proba- 
bilities are then used to obtain D(q) and the histogram 
of the probability distribution. D(q) is obtained by ap- 
plying a linear fit to logZ q in Eq. ([I]) versus log/, where 
I is the box length. The fit was performed for a range 
of I over which the function was linear. The histogram 
of the probability distribution is measured by using ex- 



FIG. 4: Exponent of the power-law fit to the probability dis- 
tribution vs. probability for the complete hull of a percolation 
cluster for several different widths. The inset shows a simi- 
lar plot for the external hull with the associated theoretical 
prediction (black line) Color online 



ponentially smaller bin sizes, e.g., the first box has size 
1/2, the next 1/4, then 1/8, etc. The exponent of the 
power-law is fit at different probabilities using 5 points 
which roughly span an order-of-magnitude in probability. 

Simulations of percolation and Ising clusters were per- 
formed for a number of system widths. Our results are 
for w = 400, 1000, 2000, and 4000, for which data was 
recorded for all systems. Small systems, w — 400, have 
w 5 • 10 5 hull sites in the cluster and large systems, 
w = 4000, have ss 5 • 10 6 hull sites. D(q) and the slope 
of the power-law fit to the probability distribution were 
obtained for the complete and external hulls of both per- 
colation and the Ising clusters. 

Fig. [3] shows a comparison between the results of the 
complete and external hulls of percolation clusters with 
the theory for the external hulls, Eq. Q. There is good 
agreement among all three for large q, which is not sur- 
prising as the complete hull fjords contribute negligibly 
to D(q). For small q there is significant disagreement 
between the complete hull and the theory as the two 
must approach different values for D(0). Previous simu- 
lations [S] have shown that D(0) increases with increasing 
width; however we see a peak at a width of 100; see Fig. 
[3] This is because there is a non-negligible fraction of 
the hull sites with probabilities less than 10 -300 for large 
widths. We expect for very large systems, if we are able 
to record all probabilities, that the complete hull D(q) 
will be nearly identical to the theory for q > 0, at q = 
there will be a jump to D(0) = 7/4, and for q < 0, 
D(q) will be ill-defined (unbounded). In comparison, for 
the external hull we see excellent agreement between the 
data and the theory Q over the entire range of D(q), 
especially for largest system sizes. 

The exponent of the power-law fit to the probability 
distribution for the complete and external hull (inset) 
are shown in Fig. [4] The exponent for the complete hull 
is —0.996 ± 0.01 over, incredibly, 50 orders of magni- 
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FIG. 5: D(q) vs. q for the complete and external hulls of 
Ising clusters for four different widths compared to theory 
from Eq. Q (black line). Inset: small q behavior of D{q). 
Color online. 

tude! We presume that the exact value of the exponent 
is — 1, which implies that D(q) is undefined for q < 0. 
Previous simulations pQ were unable to obtain this re- 
sult because the smallest probability that could be mea- 
sured, ss 1CP 10 , is still in the transient regime. The initial 
overshoot of the power for small systems corresponds to 
the probability distribution for the external hull being 
picked up by the complete hull. The power-law expo- 
nent is also obtained for the external hull, —0.93 ± 0.05, 
which is consistent with the theoretical prediction of 
-23/24 - -0.958. 

Similar results were obtained for the Ising model. 
Fig. [5] shows the comparison between the complete and 
external hulls of Ising clusters with the theory [5] for 
D(q). As with percolation, there is good agreement with 



theory for large q for both the complete and external hulls 
but significant disagreement at small q for the complete 
hull, where the Eq. Q does not apply. The probability 
power-law exponents for the complete and external hull 
are, -0.997±0.012 and -0.920±0.048, respectively for the 
Ising model. The complete hull exponent again points to 
q = as the discontinuity point for D(q); the external 
hull exponent agrees roughly with theory. 

In summary, we have described a method to obtain 
harmonic measures representing events of extremely low 
probability, which allowed us to probe the internal struc- 
ture of percolation and Ising model cluster hulls, where 
we observe the asymptotic behavior D(q) — g -1 . We are 
not aware of any theoretical prediction or explanation of 
this result. 

In future work |15j . we plan to apply the continuous 
version of this algorithm to obtain the harmonic measure 
for Diffusion Limited Aggregation (DLA) [16] for which 
there are no exact results, though there are several con- 
jectures for the form of D(q) for small and negative q 
[17. For the theory of DLA, the harmonic measure plays 
a central role, because it represents the growth proba- 
bility at every point on the cluster at a given time. The 
best current results for D(q) use iterative conformal maps 
[L8] [19], and are restricted to rather small clusters of 
w 10 4 sites. Our method will allow us to go to much 
larger sizes, w 10 7 sites. This is important because the 
slow crossover of some length scales in DLA [2D] suggests 
large clusters are necessary to approximate the scaling 
limit. Our method could shed light on the internal struc- 
ture of DLA for which little is known. 
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National Science Foundation through DMS-0553487. 
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